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APPENDIX A 


I. Introduction 


The primary goal of the research conducted under this 
grant has been and will continue the design and implementation 
of hardware and software for real-time computer graphic displays 
for cockpits. The main emphasis of the past six month period 
has been the development, simulation and testing of an algorithm 
for anti-aliasing vector drawings. 

II. Anti-aliasing of Vector Drawings 

Of great interest to the users of raster graphic display 
devices, in the removal of the adverse effects of spatial samp- 
ling. The pseudo-anti-aliasing line drawing algorithm we propose 
is an extension to Bresenham's algorithm for computer control 
of a digital plotter (I ] . While retaining the salient features 
of the original algorithm, the new algorithm does not reproduce 
a line as a sequence of disjoint line segments. The new algo- 
rithm produces a series of overlapping line segments where the 
display intensity shifts from one segment to the other in this 
overlap (transition region). True anti-aliased lines can be 
considered as having an overlapping behavior as well, but in 
these lines the rate of intensity shift and therefore the length 
of the overlap is a function of the slope of the true line. In 
this algorithm the length of the overlap and the intensity shift 
are essentially constants because the purpose of the transition 
region is an aid to the eye in integrating the segments into a 
single smooth line. 

The anti-aliasing algorithm retains the following important 
features of Bresenham's algorithm: 
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Implementation of a Simple Anti-aliasing Algorithm 


This is a brief description of an implementation of our new 
anti-aliasing line plotting algorithm for a 512x512 raster display. 
The purpose of this document is to explain how the original Bresen- 
ham algorithm was modified for anti-aliasing. 

Because the line plotting routine places pixel codes in a 
frame buffer, the intensity information which will be used in 
the discussion of the algorithm will refer to the two low order 
bits of each pixel. Full intensity, white, pixels will have 11 
as their low order bits. Black pixels will have 00 as their low 
order bits, and the algorithm calls for two intensities which 
represent two steps from black to white. The brighter of these 
two is the intermediate intensity (referred to as 66% in the pro- 
gram comments) and has 10 as its low order bits. The last 
intensity is the minimum intensity (33%)' which has 01 as its low 
order bits. This numbering scheme allows the high order bits 
to represent color and allows new pixels to be ORed into memory. 

The possibility of accidently converting an inteimediate pixel 
to a full intensity pixel by this ORing process ana the visual 
effect this causes is so slight, that the use of additional 
bits or additional code to prevent this should not be considered. 

Three program variables, FULL, IMED and IMIN, set prior to 
time generation, usually represent the full, intermediate (66%) 
and minimum (33%) intersities, respectively. IMIN will in one 
case (covered later) be set to the intermediate (66%) value; 
this is the only exception. 

There are three constants used in the algorithm; these 
should be powers of two since they are used in multiplies and 


should be implemented as shifts. The first constant is the 
number of pixels in the overlaps of the axial time segments 
(stairsteps) generated by the Bresenham algorithm. This lap 
constant (denoted LAPCON (1) in the program listings) is used 
for lines which are constructed primarily of axial moves. These 
axial lines form an angle less than tan ^(.5) with a vector 
aligned with the Ml move and are will execute an Ml move first 
(they have a negative initial decision variable) . The second 
constant, also a lap constant (lAPCON (2)), is used to set a 
long overlap used on axial lines which form very shallow angles 
with the Ml axis. These lines have long runs of Ml moves and 
the longer overlap enhances their appearance, giving a better 
approximation of the true line. The lines using the second lap 
constant have a large difference in their Aa and Ab values 
(see Bresenham) . The last constant (denote RATIO, and referred 
to as the aspect ratio) determines how great the difference 
between Aa and Ab must be to use the second lap constant. The 
following values are normally used for these constants: 

Lap constant one - 4 (pixels) 

Lap constant two - 16 (pixels) 

Aspect ratio = 32 

If these constants are changed, their relative magnitudes must 
remain the same, that is, RATIO must be the largest and lap 
constant one must be the smallest and at least equal to two. 

Execution of the algorithm begins with the normal compu- 
tations used in Bresenh2un's algorithm. The octant is established, 
the Ml and M2 moves are set, ^ and Ab are set, and the Initial 
value of the decision variable, (referred to as delta) , is 
computed . 


The next coir.putatiorb are set two test variables used to 
p>osition the overlaps r i-’et the actual intensities to be used, 
and in some cases change the value of Bresenhaun's delta. The 
two test variables (denoted ANTI2 and ANTIl) are used to 
locate the transition region's (overlap's) starting and midpoints, 
respectively. When the algorithm initially enters the transition 
region, it produces the overlap by outputting a minimum intensity 
in the M2 direction; it makes the Ml move; and puts out an inter- 
mediate intensity. At the midpoint, an intermediate intensity 
is used in the M2 direction and a minimum intensity in the Ml 
direction. Because the transition region is divided into two 
equal parts, after the computation of the midpoint test value 
(ANTIl) , the start point test value (ANTI2) is always set equal 
to twice the midpoint value (a shift left of one) . 

After the full (FULL) and intermediate (IMED) values are 
set using their respective parcuneters, the initial decision 
variable, delta (DELTA) is check for a non-negative value. 

If delta is greater than or equal to zero, then the line is of a 
diagonal type (it will contain only singular Ml moves, if any) . 
For these types of lines, ANTIl is set to its minimum value, 

-2a b, and the minimum intensity is set to the value of the 
intermediate parameter (66%). Processing then proceeds to the 
setting of ANTI2 and the generation of pixels. 

For axial type lines, a series of three cases are checked 
to set ANTIl. First Aa is checked to see if it is greater than 
or equal to the aspect ratio (RATIO) times Ab. If it is, then 
ANTIl is set to minus the long lap constant (LAPC0N(2)) times 
Ab. If this first test fails, then ANTIl is set to minus the 
first lap constant (LAPCON(l)) time Ab; this value of ANTIl is 


ANTIl l3 now compared with the initial delta computed by the 
Bresenheun algorithm. If ANTIl is less than delta, it is set to 
its minimum value, -2Ab. Now that ANTIl is set, the minimum 
(IMIN) intensity is set to the value of its parameter (33%). 

Next, ANTIl is added to delta (DELTA) to shift the laps in 
axial type lines for symmetry. After ANTI2 is set, the gener- 
ation of pixels can begin. 

The initial value of delta is compared to 7VNTI2, if it is 
less than ANTI2 then the first pixel of the line is output at 
full intensity. Otherwise, the pixel is output at the inter- 
mediate level. 

A count is now initialized using the Aa value and is decre- 
mented each time Ml or M2 loop is executed. As indicated by 
Bresenham, this value (Aa) is the number of moves necessary to 
generate the line. When the count reaches zero, line generation 
is complete. Comparing the current position with the true line 
endpoint will not always work. Although the endpoint will be 
output by the algorithm, it may be in a lap and never actually 
coincide with the position datum. 

Line generation begins now with the same loop (Ml or M2) 
that it would for Bresenham* s algorithm and will proceed until 
the count (mentioned above) reaches zero. The M2 loop is exactly 
the same as it is for the Bresenham algorithm with pixels output 
at full intensity. The Ml loop contains the additional code 
for anti-aliasing. 

When the Ml loop is entered, delta is compared with ANTI2. 

If delta is less than ANTI2, then the Ml loop performs exactly 
as specified by Bresenham and output a new pixel at full 
intensity. Since delta is usually less, the test instruction 


is generally the only new instruction executed. The FORTRAN 
listing of the algorithm uses a slightly different Ml loop to 
keep track of position# but the result Is the same. 

If delta Is greater than or equal to ANTI2# It is compared 
with ANTll. If It Is less than ANTIl, then the following 
sequence takes place. First a pixel of minimum intensity is out- 
put using the current (not updated) position plus the M2 move 
as its location. Then the position is updated with the Ml 
move and usual pixel is output but with an intermediate intensity. 
The loop completes execution by updating delta as specified by 
Bresenh£un. If delta is greater than or equal to ANTll# then 
exactly the same sequence of operations takes place# except 
that the first pixel output has an intermediate intensity and 
the second has a minimum intensity. 


Reference . 


[Ij J.E. Bresenham, "Algorithm for Computer Control of a 
Digital Plotter IBM System J. 4» 1965. 
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SUBROUTINE DRAWO( STARTX, ST ARTY. ENDX, ENDY ) 

► «*»#***#*♦♦•• •*«**•••••*••««*««•*«***»««*•«»«««*«« 




LINE PLOTTING ROUTINE USING STANDARD BRESENHAH ALGORITHM 

- THIS SUBROUTINE DOES NOT DO ANY ANTIALIASING AND 
ANY COMMENTS PERTAINING TO THAT MAY BE IGNORED. IT 
GENERATES ONLY FULL INTENSITY PIXELS (INTENS(l)) 


CALL PARAMETERS 


STARTX. STARTY » X 
Et4DX.E14DV a X 


Y COORDINATES OF LINE START POINT 

Y CORRDINATES OF LINE END POINT 


»*•>**'««**«*««»•«'*«««*•**«#••«*«»««**«**«#*««»«****«*«*•*«**•*«'*« 

IMPLICIT INTEGER*2(A-2) 

COMMON RATIO. LAPC0N(2). INTENS ( 3 ) . FRAME ( 51 2, 512). DlAC 

*«****•*••*«*«"*««*****•«*»*««*«««•»*«««»»««»«*****•***«»•««.*«*•« 

COMMON BLOCK PARAMETERS 


FRAhF 

= PICTURE ARRA> 

1 NTENS 

= INTENSITY TABLE (VALUES AS FOLLOWS) 

1 FULL INTENSITY PIXEL CODE 

2 = 66% INTENSITY PIXEL CODE 

3 » 33% INTENSITY PIXEL CODE 

RATIO 

= ASPECT RATIO FOR ANTIALIASING ROUTINE. USED TO 
DETERMINE SHALLOW AND STEEP (NEAR 4b DEG. ) LINES 
WHICH USE A LONGER LAP (SEE LAPCON) 

LAP cot; 

= PIXEL LA® CONSTANT TABLE, USED TO SET THE NUMBER 
PIXELS IN THE TRANSITION (LAP) REGION (VALUES AS 
FOLLOWS ) 

1 = STANDARD LAP (INTERMEDIATE SLOPE LINES) 


2 a long lap (FOR SHALLOW AND STEEP LINES) 

COMPUTE 

DELTA X AND Y AND SETUP OCTANT 


DELX » ENDX - STARTX 

DFLV = ENDv - STAR TV 

DELXV = IABS<DELX> - IABS(DELV ) 
OZTAtJT 1 OR 2 

IF (DELX CE 0 AND DELY CE 0) THEN 
M2X = 1 

M2Y « 1 


OCTANT 1 


IF (DELXY CE 0) THEN 
DELA = DELX 
DELB a DELY 


./RAWO 


0013 


MIX » 1 




0014 


MIY • 0 





C 

OCTANT 2 




0015 


ELSE 




0016 


DELA « DELY 




0017 


DELB *> DELX 




0018 


MIX ■ 0 




0019 


MIY » 1 




0020 


END IF 





C 

octant 3 or 4 




0021 


ELSE IF (DELX LT. 0 AND DELY 

GE 

0) 

THEN 

0022 


M2X = -1 




0023 


M2Y ■ 1 





C 

0-CTA.JT 4 




0024 


IF (DELXY. CE. 0) THEN 




0025 


DELA ' -DELX 




0026 


DELB «= DELY 




0027 


MIX » -1 




0028 


MIY « 0 





C 

o:tant 3 




0029 


ELSE 




0030 


DELA » DELY 




0031 


DELB “ -DELX 




0032 


MIX * C 




0033 


MIY * 1 




0034 


END IF 





c 

0-:TANT 5 OK e 




0035 


ELSE IF (DELX LT 0 AND DELY 

LT 

0) 

THEN 

0036 


M2X -1 




0037 


M2T -1 





c 

OCTANT 5 




0038 


IF (DELXY CE 0) THEN 




0039 


DELA = -DELX 




0040 


DELB » -DELS 




0041 


MIX » -1 




0042 


MIY * 0 





c 

o*:tant 6 




0043 


ELSE 




0044 


DELA » -DELY 




0045 


DELD « -DELX 




004t 


MIX 0 




0047 


MIV » -1 




0048 


END IF 





c 

o:tant 7 OR 8 




0049 


ELSE 




0050 


M2X » 1 




0051 


M2Y * - . 





c 

OCTANT e 




0052 


IF (DELXY CE 0) THEN 




0053 


DELA DELX 




0054 


DElB = -DELY 




0055 


Mix = 1 




0056 


MIV » 0 





c 

OrTANl 7 




0057 


ELSE 




0058 


DELA «= -DELY 




0059 


DELD » DELX 





DhAWO 


0060 

0061 

0062 

0063 


0064 

0065 

0066 

0067 

0068 
0069 


0070 

0071 

0072 
0072 

0074 

0075 


0076 

0077 
0076 

0079 

0080 
0081 
0082 

0083 

0084 

0085 

0086 
0087 
0086 
008*5 

0090 

0091 


C 

C 

C 

C 

C 


END IF 


END IF 


MIX 

MIY 


■ 0 
- -1 


SETUP LINE DRAWING ALGORITHM 


SET UP CONSTANTS 

DL 20 » 2*DELB 

DEL2A0 * 2«(DELD - DELA) 

DELTA « DEL2D - DELA 

OLDX * STARTX 

OLDY * STARTY 

FULL * INTENS(l) 


C** ♦■•♦**♦■»*♦•»••♦•»♦♦♦■•••♦♦**♦♦♦•♦#♦*♦•♦■*•♦*•♦••****♦****•♦* **************** i 

c 

C DIAGNOSTIC ROUTINE 

C 

c ♦•*****•♦♦♦♦*«***«*♦«♦**♦••♦♦•*♦♦♦ ♦**♦♦•***«♦***♦***♦**♦'♦•*♦*♦**♦*♦***♦♦* 

IF (DIAC CT 0) THEN ' IN DIAGNOSTIC MODE 

WRITE <6, 2030) 

203C FORMAT (' DRAWO SUBROUTINE - STANDARD BRESENHAM ' ) 

WRITE<6. 2032) DELA. DELB. DELTA 
2032 FORMAT< ' A*'. 14. ' B='. 14. ' DELTA*'. 15) 

END IF 


C 


c:«vr*«**«*«««««*«**««-»****«*««*«**«*»«««««««««««*«««***«*******«***** 

c 

C DRAW THE LINE 

C 

(:«**«*««««*«*•«*«*•«*««*««*«*«•**«*««••««*«*•««•«•*««•«•**••«•••«**« 

C OUTPUT THE STARTING POINT 

FRAME (OLDX. OLDY) » FULL 
C DIAW ,,'E REMAINDER OF THE LINE 

100 IF (DELA CT 0) THEN • DELA » NO OF POSITIONS IN 

IF (DELTA LT 0) THEN 

C Ml MOVE 

OLDX ■ OLDX ♦MIX 

OLDY « OLDY ♦ MIV 

FRAME (OLDX. OLDY) « FULL 
DELTA - DELTA ♦ DEL2B 

ELSE 

C M2 MOVE 

OLDX = OLDX ♦ M2X 

OLDY « OLDY ♦ M2V 

FRAME (OLDX. OLDY) « FULL 
DELTA » DELTA ♦ DEl2AB 

END IF 

del A - DELA - 1 • DECREMENT POSITION COUNT 

GOTO 100 

END IF 


LIh 


DR AMO 



C 



c 

LINE DRAWING COMPLETED 


c 


0092 


RETURN 

0093 


END 


000 1 


0001 ' 
00 Cl' 


0004 

0005 

0006 


SUDROOTINE DRAWl (STARTX. STARTY. ENDX. ENDY) 

C **«««««»«***»*«*«•*««««««««*«««*•««««»««•««•• »««««■••«••**«»«««««»«««•» 

c 

C LINE PLOT! INS ROUTINE USING ANTIALIASING 3RESENHAM ALGORITHM 

C DEVELOPED BY E J DUNNING THIS ROUTINE USES THE SIMPLE 

C VERSION OP THE ANTIALIASING ALGORITHM WHICH ONLY CREATES 

C TRANSITION REGION LAPPING ON LINES WHICH ARE PRIMARILY 

C AXIAL IN NATURE THESE LINES ARE IDENTIFIED BY A NEGATIVE 

C INITIAL DELTA VALUE COMMENTS IN THIS ROUTINE PERTAINING 

C TO ANTIALIASING STEEP LINES MAY BE IGNORED. THE APPEARENCE 

C OF THE STEEP LINES IS ENHANCED. HOWEVER. BY THIS ROUTINE BY 

C USING 66*/. INTENSITY PIXEL VALUES IN THE MI TRANSITIONS 

C 

c BY E JACK DUNNING 

C 

C *«»«’*«•*****«****«««'***«««'• ****** ««'«•«*«»«»«««••*«••«««»**•«*«•**•*«'»» 

c 

C CALL PARAMETERS 


C STARTX, STARTY = X V COORDINATES OF LINE START POINT 

C ENDX.ENDV = X S, Y CORRDINAIES OF LINE END POINT 

C 

C*'*f**«-*««'»«*'*»******-*'**#'***** ♦*•♦♦♦♦♦♦♦**♦**♦♦**** *♦**♦**♦♦*♦*■«>* ' i « ' V . 

IMPLICIT INTEGER^OtA-Z ) 

COMMON RATIO, LAPC0N(2). INTENSiO) , FRAME(5l2. 512>. DIAC 

♦♦• ******♦♦<!*♦ »*♦*♦**#*•***♦****#**♦♦••*♦# 

c 

C COMMON BLOCK PARAMETERS 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


FRAKE = PICTURE ARRAY 

INTENS - INTENS.ITY TABLE (VALUES AS FOLLOWS) 

1 = FULL INTENSITY PIXEL CODE 

2 66*/. INTENSITY PIXEL CODE 

3 = 33 % INTENSITY PIXEL CODE 

RATIO = ASPECT RATIO FOR ANTIALIASING ROUTINE. USED TO 

DETERMINE SHALLOW AND STEEP (NEAR 45 DEG ) LINES 
WHICH USE A LONGER LAP (SEE LAPCON) 

LAPCON = PIXEL LAP CONSTANT TABLE. USED TO SET THE NUMBER 
PIXELS IN THE TRANSITION (LAP) REGION (VALUES ' 
FOLLOWS ) 

1 * standard LAP (INTERMEDIATE SLOPE LINES) 

2 ■= LONG LAP (FOR SHALLOW AND STEEP LINES) 




c 

C COMPUTE DELTA X AND V AND SETUP OCTANT 

C 

c «••«*«•*«•*•«««»«*«••«*••«•«•*«••«**•*« ••««•* 

DELX = ENDX - STARTX 
DELY = ENDY - STARTV 
DELXY = lABS(DELX) - lABS(DELY) 

C O'- T ANT 1 OR r 



DRAW! 


0007 

0008 

0009 

0010 
001 1 
0012 

0013 

0014 

001 
00 1 o 

0017 

0018 

0019 

0020 


0021 

0022 

0023 


0024 

002 

0026 

GOO” 

002F 

002 ^ 

0030 

0031 

0032 

0033 

0034 


003 ^ 

0036 

0037 


0038 

0039 

0040 

0041 

0042 

0043 

0044 
004 5 

0046 

0047 
0046 

0049 

0050 

0051 

0052 

0053 


IF (DELX GE 0 AND DELY OE 0) THEN 
M2X = 1 
M2Y « 1 

C OCTANT I 

IF (DELXY CE 0) THEN 
DEL A - DELX 
DELB = DELY 
MIX ■ 1 
MIY « 0 

C OCTANT 2 

ELSE 

DELA - DELY 
DELD = DELX 
MIX - 0 
MIY * 1 

END IF 

C OCTANT 3 OR 4 

ELSE IF (DELX LT 0 AND DELY GE 0) THEN 
M2X = -1 
M2Y = 1 

C OCTANT 4 

IF (DELXY GE 0> THEN 
DELA *= -DELX 
DELP »= DELY 
MIX = -1 
MIY = C 

C OCTANT 3 

ELBE 

DELA = DELY 
DELD = -DELX 
MIX = 0 
MIY » 1 

END IF 

C OCTANT 5 OR 6 

ELSE IF (DELX LT 0 AND DELY LT O) THEN 
M2X = -1 
M2Y = -1 

C OCTANT 5 

IF (DELXY GE 0> THEN 
DELA -DELX 
DELB -DELY 
MIX = -1 
MIY = 0 

C OCTANT 6 

ELSE 

DELA = -DELY 
DELD = -DELX 
MIX « 0 
MIY = -1 

END IF 

G OCTANT 7 OR 0 
ELSE 

M2X » 1 
M2Y - -1 

C OCTANT 8 

IF (DELXY GE 0) THEN 
DELA « DELX 



DRAUl 


0054 

0055 

0056 

0057 
005B 

0059 

0060 
0061 
0062 
0063 


0064 

0065 

0066 
0067 
0066 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 
007B 

0079 

0080 
0081 

0082 

0083 

0084 

0085 

008b 

0087 


0086 

0089 

0090 

0091 

0092 

0093 


C 


OCTANT 7 


ELSE 


END IF 

END IF 


DELB “ -DELY 
MIX » 1 
MIY - 0 


DELA » -DELY 
DELB « DELX 
MIX « 0 
MIY - -1 


c 

C SET-UP LINE DRAWING ALGORITHM 

C 

C ***«•***«»'»«•»*««««•«»«'•««***** ••«*«*'»•»*****•***•»»*«««'*«**«»*****'•'*< 

C SET-UP CONSTANTS 

DEL2B B 2«DELB 

DEL2AB = 2* (DELB - DELA) 

DELTA = DEL2B - DELA 
OLDX = STARTX 
OLDY = STARTY 
C SET-UP ANTIALIASING 

FULL = INTENS(l) 

IMED = INTENS(2) 

IF (DELTA GEO) THEN ' DIAGONAL TYPE LINES 

ANTIl = -DEL2B ! DEFAULT VALUE 

IMIN * INTENS(2) ! SET MIN TO MED INTENSITY 

TYPE = 3 


ELSE 'AXIAL TYPE LINES 

IF ( (RATIO*DELB ) . LE, DELA) THEN ! SHALLOW AXIAL 

ANTIl = -LAPC0N(2> * DELB ! LONG LAP 

TYPE = 1 


ELSE ' STANDARD AXIAL 

ANTIl = -LAPCON(l) * DELB ! STANDARD LAP 

IF (DELTA GE ANTIl ) ANTIl = -DEL2B 

' MINIMUM TRANSITION REGION (LAP) LENG' 


TYPE B 2 


END IF 

IMIN = INTENSO) ' STANDARD MINIMUM INTENSITY 

DELTA = DELTA ♦ ANTIl ' CORRECT SYMMETRY 


END IF 

ANTI2 = ANTIl * 2 




n 

L. 

C 

c 

c 


2030 

2032 


DIAGNOSTIC ROUTINE 


IF (DIAG GT 0) THEN ' IN DIAGNOSTIC MODE 

WRITE (6. 2030) 

FORMAT(' DRAWl SUBROUTINE - ANTIALIAS ALG NO 1') 
WRITE (6. 2032) DELA. DELB. DELTA, ANTIl. ANTI2. IMED, IMIN, TYI 
FORMAT(' a='.I4,' Bb'.I4, ' DELTA* ' . 1 5. ' A1-',I5, 

' A2*',I5, ' MED*', 14, ' MIN*', 14, ' TYPE*', 12) 

END IF 


DRAUl 


C« ««««««*«««*««««»*«*««»*«**«««*«»« ««««««»««»«««-*«««««»«*»««««»»«««»»««! 


0094 

0095 

0096 

0097 

0098 

0099 

0100 


0101 

0102 

0103 

0104 

0105 
OlOfc 


0107 

0106 

0109 


0110 

01 1 1 
0112 
0112 

0114 

0115 

0116 
Oil” 
0116 
01 1<= 
0120 
0121 
0122 
0123 


0124 

0125 


1 

C DRAU THE LINE I 

C 

C***^************ «■«««***«*»**«**««***«««**«****«*«»««»»*»*«*•**««**«««* 

C OUTPUT THE STARTING POINT 

IF (DELTA. LT. ANTI2) THEN 

FRAKE(OLDX,OLDY) * FULL 

ELSE 

FRAME ( OLDX. OLDY) = lOR ( FRAME ( OLDX, OLDY ) . IMED ) 

END IF 

C DRAW THE REMAINDER OF THE LINE 

100 IF (DELA CT. 0) THEN • DELA = NO. OF POSITIONS IN LI 

IF (DELTA. LT 0) THEN 

C Ml MOVE 

NEWX = OLDX + MIX 
NEWY = OLDY + MIY 
IF (DELTA. LT ANTI2) THEN 

FRAME (NEWX. NEWY) = FULL 
ELSE IF (DELTA. LT. ANTI 1 ) THEN 
FRAME (NEWX, NEWY) = 

* IOR( FRAME (NEWX. NEWY). IMED) 

FRAME ( 0LDX+M2X. 0LDY+M2Y) = 

* lOR (FRAME (0LDX+M2X. 0LDY+M2Y). IMl 

ELSE 


FRAME (NEWX. NEWY) = 

* I0R(FRAME(NEWX, NEWY). IMIN) 

FRAME ( 0LDX+M2X. 0LDY+M2Y) = 

« lOR (FRAME (0LDX+M2X. 0LDY+M2Y), IME 

END IF 

OLDX NEWX 

OLDY = NEWY 

DELTA = DELTA + DEL2D 

ELSE 


C 


C 

C 

C 


M2 MOVE 

OLDX = OLDX + M2X 
OLDY *= OLDY '• 

FRAME (OLDX. OLDY) = FULL 
DELTA = DELTA DEL2AB 

END IF 

DELA = DELA - 1 • DECREMENT POSITION COUNT 

GOTO 100 

END IF 

LINE DRAWING COMPLETED 

RETURN 

END 
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MULTIPLIER ACCUMULATOR CAROS AND 
COORDINATE TRANSFORMATIONS 


Alireza F. Faryar 


Sarah A. Raj ala 


I . INTRODUCTION 


The purpose of this report is a study of four pre- 
viously built Multiplier Accumulator Cards (MAC) and their 
application to coordinate transformations. These cards are 
part of a real-time raster graphic display system used for 
generating raster graphic algorithms to be used in the cock- 
pits of aircrafts. In this case the MAC will be used to 
perform coordinate transformations* such an example is illus- 
trated in Figures 1 and 2. For real time operation this 
transformation must be very fast, and in Section II, it is 
shown that it can be done in 6.7 micro seconds. First, 
however, the hardware is discussed. 

A. Multiplier Accumulator Cards Hardware t 

These cards are powerful processing elements 
each containing a fast multiplier chip, a 32 bit ALU, input- 
output memories, and a microprogrammed controller. A simple 
block diagram is shown in Figure 3. In the following sections, 
different parts of the card and their performance are discussed, 
sed. 

1. Busing: 

Communication with each card is provided via 
three different buses. Address bus. Data bus and System 
Function bus. 

(a) Address bus: 

The address bus is a 24 bit wide and is 


used as follows: 


ABOO through AB07 provide addresses corresponding to X and Y 
when initial data is being written into the input memories. 
AB8 specifies if we are writing into X or Y memory. AB16 
controls the state of the latches (B1 to B4) . MPIX>D-L set 
low lets B7, B8« d, C8, D7 and D8 be loaded from the data 
bus. If every input to B12 is in the proper mode MPLOD-L 
will go low, when CLK goes high and if AB16 is high. AB18 
and AB19 address specific cards as follows: 


A6 

18-19 

00 

01 

10 

11 


Address Card No. 

0 

1 

2 

3 


( 


AB21-AB23 are an enable signal to (B15)*. 

(b) Data bus: 

The data bus is 32 bit bus used as fol- 


lows: 


DB00-DB23 are used to: 

• Load the counters B7, B8, Cl, C8, D7 and D8 with 
initial addresses of X, Y, Z when a function is 
to be performed. A new address is provided to 
these counters whenever MPLOD-L goes low (as di- 
rected in (a) . 

* If one is writing data into input memories, X and 
Y, MPLOD-L is high (latches are closed) and data 
are carried in through A1 to A8 bus drivers to the 
RAMS (this time this bus is providing input data) . 


* 


Chip number used in wiring diagram. 



If a function is performed, the 32 bit output of Z 
memory is written on the entire 32 bit data bus, 
through A1-A8. The state of bus drivers A1-A8 is 
controlled by RDIM-H provided from (B14) . When 
RDIM-H is high, the output of Z RAM is written on 
the bus. When it is low the data or the address 
on the bus is being read into the card. 

(c) System Function bus: 

The system function bus is a six bit bus 

consisting of: 

Reset-Clock-F3-F2-Fl-F0 

F0,F1,F2 initialize micro-prograun counter through address 
memory (C12) . This counter provides the address to the micro- 
progreun storage ROM's. 

When F3 is low: 

If AB8 = 0, AB16 = 0 input is being written into X RAM. 

If AB8 « 1, AB16 “ p input Y is to be written. 

If AB8 s 0 and AB16 « 1 the state is changed, RDIM-H > 0. 

The address will be provided to the address counters through 
the data bus. CLK and RESET are provided by the main system. 

2. Microprogram Counter: 

This section consists of a 256 bit bipolar PROM 
(Cl 2) and two synchronous four bit binary counters. The carry 
look ahead circuitry of the counters has made it possible to 
cascade them to get an eight bit synchronous output (C14 t C15 
FO, Fl and F2 of the function code will address the PROM. 

It's output will be an appropriate address to the micro-code 
memory for executing a specific function. The counters are 
loaded with this address and as long as MAC is in this routine 


MPINC-H (RUN-H) signal (one of the outputs of micro -prog reun 
memory) enables the counters to increment with CLK. 

3 . Microprogram Memory 

The microcode* is stored in four 2048 bit (256x8) 
bipolar PROMs in order to construct a 256x32 bit memory. 
Address to this memory is provided by the microprogram coun- 


ter 

. The output is 

used to provide appropriate signals for 

every function to be executed. 



4. 

Card Decoder 





The card 

decoder is a 

three to eight line 

decoder 

multiplexer (B15) . 

One output 

is used in each card 

to indi- 

cate if the corresponding card 

has been selected. The inputs 

to 

the card decoder 

and corresponding pin number in 

each card 

are 

shown below: 





AB 

AB 




75 71 2T 

T9 IF 

pin# 

MAC# 


Oil 

00 

7 

0 


Oil 

01 

9 

1 


Oil 

10 

10 

2 


Oil 

11 

11 

3 

For 

example, if AB 

18 and 19 are 10 then pin number 

10 of (B15) 


in card number 2 will go low. 
5. Bus Receiver-drivers 


Except for the data bus and the card busy signal^ hex- 
inverter interface elements are used to let a MAC communicate 
with the system buses. For the data bus and card busy signal 
Tri-state quad bus transceivers have been used. 


* Lists of the microcode are given in Appendix A. 


6 , Latches 

All latches are Hex D-type flip-flops with clear. 
Four of them (R1 to B4) are used to provide either starting 
addresses to Y and Z memories# when MPLOD-L is low# or to 
latch the input to the counters# when the data bus carries in- 
put data to input memories (MPLOD-L high) . 

Latches G8 to G12# E5 to E8 and part of El2 are used to 
provide the outputs of the multiplier chip and Z RAM to the 
arithmetic logic unit# whenever appropriate. They are con- 
trolled by signals provided by the microcodes. 

7. Row-counters# Column-counters and Selectors 

This circuit is explained for the X input. Similar 
circuits are used for Y and Z memories. 

(a) If data is to be written in the X input 
memory# the address will be ready on the 
address bus and data will be provided 
by the data bus. In this case# selectors 
(B5 and B6) will select the address bus 
as address to the X RAM. This happens 
when WRIM-L is low. The data on the data 
bus are now inputs to X and they are 
prohibited from appearing on inputs of 
counters (B7 and B8) by MPLOD-L signal 
being high. Instead# they are read 
directly from the outputs of Bus-Trans- 
ceivers into the X RAM. Signals XWREN-L 
specify whether the data is to be %n:itten 
in X or Y memory. 


(b) If the data la ready in Y and X RAMS and 

a function ia to be per formed « MPLOD-L will 
90 low providing a atarting addreaa to B7 
and B8 countera. WRIM~L will be high ao 
that the output of the countera will addreaa 
the X-RAM. The counter will count with 
XLXX>D*‘L aignal provided by microcode. 

(c) To read the data out and write it on the 
data bua» ROIM-L will go low and ao will 
aelect the addreaa bua aa addreaa to the Z 
RAM. Thia will provide the 32 bit long 
output of Z RAM to the bua tranaceivera A1 
to A8. Since TE > RDIM-H ia high the out- 
puta will be written on the data bua. 

8 . Memory 

There are three different typea of memory uaed in 
each card, X, Y and Z, all the memory ia made from TTL 256 X 
4 bit fully decoded random acceaa memory (93L422) . The eight 
bit addreaa to each of memory ia uaually treated aa two 4 bit 
fields of row and column addreaa. For example X(7,3) ia atored 
in location 01110011. 

(a) X-input memory 

Thia ia a 256 X 16 memory. The output ia 
alwaya enabled. The addreaa to it ia pro- 
vided either by the addreaa bua (in write 
mode when XWREN-L ia low) , or by the addreaa 
countera B7 and B8 (when executing a function, 
XWREN-L ia high) . 

The input to X ia provided by the data bua, 
when it ia in read mode. 

In thia caae the following aignala are pro- 
vided : 

RDIM-H low 

MPIX)D-L high B1 to b 4 are latched 
ABB - 0 therefore XWREN-L - low 


The output of this memory Is connected to 
multiplier chip. 

(b) y - Input memory 

This 256 X 16 RAM is addressed similar to 
X input memory. The only difference is that 
the output enable signal is not held at a 
fixed level. The multiplier chip MPY>'16AJ 
(Fl) receives Y inputs from the same pins 
where it delivers half of its output. The.re- 
fore# a signal called Tril is used to specify 
if the chip is reading or writing. The in- 
verted Tril-H is used to enable the output of 
Y memory which is directly connected to multi- 
plier chip. The input to Y is provided by 
the data bus as it was to X. 

(c) Z - output memory 

The Z output monory is a 256 X 32 bit random 
access memory. The address is provided either 
by the address counters (D7 ft D*) when it is 
in reading mode# or by address bus when it 
is in write mode. It is in write mode when 
the write enable signal, ZWREN-L, is low. 

Its input is the output of arithmetic logic 
unit. The output enable signal is grounded 
so the output is always enabled if ZWREN-L 
is high. The output is provided both to the 
bus tranceivers (A1-A8) and the latches (E5- 
E8) . The latter connects Z output back to 
the ALU. 

9. Multiplier chip MPY-16J 

The multiplier chip is a sinole chip on each MAC. 
This chip simply multiplies t%fO 16 bit binary numbers and pro- 
vides a 32 bit long product at its output. This is done in less 
than 200 n-sec. 

In order to save space on the chip, one input shares 


pins with 16 bits of the output (LSP-out) . A signal (TRIL) 


controls the flow of the data on these 16 bits. When Ttil is 
high the output of the Y-RAM i« provided to the multiplier and 
the chip will read in the data. When Tril is low the least 
significant protion of output will be ready at the output. 

The second input is provided by the X RAM. Every time a multi- 
plication is done the output will go to the arithmetic logic 
unit for further arithmetic operatiots. Since the least signi- 
ficant portion of output and the Y input use the same pins, 
the output may not be directly connected to the arithmetic 
logic unit. The connection is made through a series of latches. 
Both multiplier chip and its latches are clocked by (MULAEN- 
H & MULCLK-H) provided in the microcode. For specifications 
see Appendix C. 

10. Arithmetic Logic Unit 

The ALU consists of eight ALU chips and three look- 
ahead carry generators (74sl01 6 74sl82) . These tv>o chips 
together can perform high speed arithmetic operations. Al- 
together^ it is a 64 bit input-32 bit output ALU with a full 
carry look-ahead scheme. The inputs are provided by Z-RAM 
and multiplier chip. This circuit performs 16 binary arithmetic 
operations, with the functions and modes of operation selected 
through 5 inputs (one for mode of operation, 4 for function 
selection), which are provided in microcodes (ALUFO, 1,2,3, 
and ALUMO and ALUCO) . The output is directly connectctd to 
the Z-RAM. 

For ease of understanding a flow-chart of read and 
write operations is given in figure 4 and a timing diagram is 
shown in figure 5. 


1 . Subroutines 


There are four different subroutines considered In 
a MAC. Every subroutine Is called by an appropriate function 
code. The MAC halts after each operation and Is ready for a 
new command. In the following routines X(0,0)» Y(0»0), and 
Z(0,0) are all offset by the starting address. 

2. Vector dot product 

The following eight dot products are formed In 6.7 
microseconds . 

Z(0,0) = X(0,0)*Y(0,0)+X(0,1)*Y(1,0)+X(0,2)*Y(2,0)+X(0,3)*Y(3,0) 
Z(1,0) » X(1,0)*Y(0,0)-I-X(1,1)*Y(1,0)+X(1,2)*Y(2,0)+X(1,3)*Y(3,0) 

0 

Z(7,0) = X(7,0)*Y(0,0)+X(7,1)*Y(1,0)+X(7,2)*Y(2,0)+X(7,3)*Y(3,0) 

3. Perspective multiplication 

The following sixteen pairs of products are formed in 
6.7 microseconds. 

Z(0,0) = X(0,0)*Y(0,0) 

Z(0,1) = X(0,1)*Y(0,0) 

Z(1,0) « X(1,0)*Y(1,0) 

Z(l,l) » X(1,1)*Y(1,0) 


Z(15,0) = X(15,0)*Y(15,0) 
Z(15,l) = X(15,1)*Y(15,0) 


f 


4. Weighted Sum 


The two following 4X4 weighted sums are formed in 

6.7 microseconds. 

2(0,0) - I I X(i,j)*Y(i,j) 
i-0 j»0 

2(1,0) -IE X(i,j)*Y(i-4, j) 
i-0 j-0 

5. Vector transformation 

A 1X4 transformed matrix is formed by multiplying a 
4X4 matrix (transformation matrix) with the 1X4 original vector, 
in 3.5 microseconds. 

2(0,0) = X(0,0)*Y(0,0)+X(0,1)*Y(1,0)+X(0,2)*Y(2,0)+X(0,3)*Y(3,0) 
2(0,1) » X(0,0)*Y(0,1)+X(0,1)*Y(1,1)+X(0,2)*Y(2,1)+X(0,3)*Y(3,1) 

2(0,3) - X(0,0)*Y(0,3)+X(0,1)*Y(1,3)+X(0,2)*Y(2,3)+X(0,3)*Y(3,3) 

6. Applications 

As was previously mentioned one special application of 
these cards is the real time transformation from figure 1 to 
figure 2. This is done by transforming starting and ending 
points of each line in figure 1, using vector transformation 
routine. Prior to multiplication, a transformation matrix is 
created by concatenating a rotation matrix (figure 6. a) and 
a translation matrix (figure 6.b) . This multiplication will 
transfer a point in three dimensional space to another point 
regarding new location and position of the aircraft. A separate 
routine is used to create a two-dimensional representation of 
the scene. Finally, the picture is drawn using a line drawing 


routine. 


Another example » one can use a weighted summing to perform 
convolution. This can be done by proper selection of X(0#0) , 
Y(0,0), 2(0,0). 



IV 
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Figure 1 


Exanple of the displayed picture by the system on board. 
In this case the aircraft Is In proper path for landing. 
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Appendix B.l 
Microcode Table 
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Multiplier Chip Specification 
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Figure 3. Output Delay Versus Temperature 


Figure S Multiply Time Versus Temperaturi 
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Figure 4. Three State Delay Versus Temperature 
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Figure 1. 4*Bii Two't Compitmtnt Range 

Two*i compitmtnt nctation it tsptcially uttiul to many computti tvtttmt. Ii ofitrt trit advantage of havirtg only a tingle reprt 
temaiion for the number aero, at opposed to two for tign magnitude and orte*t complement tyttemt AdditiortaMy, tt preduoei 
the ute of **tubtractort** - potititie or negative numbers may be added to orte another without any regard for the sign of the 
numbers: the result will always be correct in two*s complement rtoiation. 

Although positive number representation is the same for both sign-magnitude and two's complement, rtegatwe numbers are 
determined by the followirtg equation: 


(Two's Complement! 


2*^ - |X| 


where |X| is the magnitude of the desired negative number and N is the total number of bits used in the two's complement field. 
iiKluding the sign bit. for integers (N ■ 1 for fractions!. 

Although the preceding might seem to be a cumbersome funaion one must perform before every negation, it turns out that it is 
easily implemented with hardware. The equivalent of the above equation in hardware is simply the inversion of all bits of D(l« 
including the sign bit. plus the addition of birsary *1* to the LSB. The same procedure reapplied to the two's complement nsimbef 
yields IXl. Inversion is quite straightforward and the addition can generally be performed with the typically unused Carry-In input 
in the LSB adder. 

FRACTIONAL/INTEGER MULTIPLICATION 

The MFY-Series multipliers may be used m either a fractional or integer mode • the difference is conceptual. For eaample. uung 
a 4-bit case, the multiplier does not know lor care! whether it is performing the multiplication 6 a (•2!« *12 or 16/B! a l-2/Bl 
• -12/B4: the input and output bmary fields will be the same. Fractional multiplication (usmg fields as defined m the previous 
specificationsl offers the advamage of more c o nvenient single precision usage. The MSB is that which is dosesi to the birsary 
point in fractional representation (the LSB for integer representation!. The fraaiorsai notation is additionally the most convenieni 
when implemeniing a floating point multiplicaiion system 



▼nw ntaenvea tne hiomt to CMauoe anoouers and sac civic at lona without noticc this 
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